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We present the results of charged particle fluctuations measurements in Au + Au collisions at 
y/SNN = 130 GeV using the STAR detector. Dynamical fluctuations measurements are presented 
for inclusive charged particle multiplicities as well as for identified charged pions, kaons, and protons. 
The net charge dynamical fluctuations are found to be large and negative providing clear evidence 
that positive and negative charged particle production is correlated within the pseudorapidity range 
investigated. Correlations are smaller than expected based on model-dependent predictions for a 
resonance gas or a quark gluon gas which undergoes fast hadronization and freeze-out. Qualitative 
agreement is found with comparable scaled p-|-p measurements and a HIJING model calculation 
based on independent particle collisions, although a small deviation from the scaling dependence 
expected from this model is observed. 

PACS numbers: PACS number: 25.75.Ld 

Keywords: Relativistic Heavy Ions, Event-by-event fluctuations. 



A key question of the heavy ion program at the Rel- whether the hot matter produced in the midst of heavy 
ativistic Heavy Ion Collider (RHIC) is to understand ion collisions undergoes a transition to and from a quark 

gluon plasma (QGP) phase before it hadronizes. One of 
the most striking signatures of such a QGP-HG (hadron 

gas) phase transition could be a strong modification in 
* URL: |wHw ■ star ■ bni ■ gov | the fluctuations of Specific observables measured on a per 
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collision basis, i.e. event by event Most often 

discussed are mean transverse momentum fluctuations 
(temperature fluctuations) and particle multiplicity fluc- 
tuations. For the latter, predictions range from enhanced 
multiplicity fluctuations connected to the production of 
QGP droplets and nucleation processes in a first order 
QGP - HG phase transition, to a strong suppression of 
fluctuations as a consequence of rapid freeze-out just af- 
ter the phase transition |^|^. In this case, final state 
values of conserved quantities, such as net electric charge, 
baryon number, and strangeness would not be strongly 
modified from their values in the QGP stage. Due to 
the large difference in the degrees of freedom in the QGP 
and HG phases, measured fluctuations, of the net elec- 
tric charge in particular, could be reduced by a factor 
ranging from 2 to 4 if a QGP is produced 0, . The 
frequency of production and size of QGP droplets may 
critically depend on the collision impact parameter. Cen- 
tral collisions are generally expected to lead to larger and 
more frequent QGP droplet production. An increase in 
the size and production frequency of QGP droplets with 
increasing collision centrality might then be signaled by 
a sudden change in the fluctuations of produced particles 
such as anti-protons and kaons fsj, as well as pious. 

In this paper, we report on a measurement of charged 
particle multiplicity fluctuations as a function of colli- 
sion centrality in Au -|- Au collisions at an energy of 
y^SNN = 130 GeV. We study event-by-event fluctua- 
tions of conserved quantities at near-zero rapidity in the 
center-of-mass rest frame (mid-rapidity). Specifically, we 
discuss fluctuations in the difference of the number of pro- 
duced positively and negatively charged particles (multi- 
licities) measured in a fixed rapidity range, defined as 



i 



(1) 



where and are multiplicities of positive and nega- 
tive particles calculated in a specific pseudorapidity, and 
transverse momentum range. The notation "(O)" de- 
notes an average of the quantity O over an ensemble 
of events. The method used to calculate the averages 
(A^+) and {N-), which vary with collision centrality, is 
described in the following (see Eqs. I5I1UI ). We consider 
fluctuations in the production of all charged particles, 
and N- (mostly pions) as well as specific cases of proton 
and anti-proton, Np and iVp, and positive and negative 
kaon, Nk+ and Nk-, fluctuations. The former amounts 
to a measurement of net electrical charge fluctuations, 
whereas the latter corresponds to measurements of net 
baryon number and net strangeness fluctuations. The 
method used to calculate this and other observables used 
in this work is described in the following. 

A difficulty inherent in the interpretation of measure- 
ments of multiplicity fluctuations is the elimination of 
effects associated with uncertainties in the collision cen- 
trality, often referred to as volume fluctuations. Event- 
by-event impact parameter variations, in particular, in- 



duce positive correlations in particle production which 
do not depend on the intrinsic dynamical properties of 
the colliding system, but rather simply reflect changes in 
the number of collision participants. Fluctuations in the 

difference of relative multiplicities, , defined in Eq. 1, 

are however free from this problem. This analysis is thus 
restricted to the study of such relative multiplicities. As 
shown in 0) ^^h can be readily translated into observ- 
ables D, and ujq, discussed by other authors ^lEBl' Its 
relation to the two-particle density is discussed below. 
We will additionally study the behavior of relative mul- 
tiplicities, i^-i , and other quantities of interest defined 

in this paper as a function of the collision centrality es- 
timated on the basis of the total charged particle multi- 
plicity measured in the pseudorapidity range \ri\ < 0.75 
in order to identify possible changes in the fluctuations 
with collision centrality. 

The magnitude of the variance, vjf , is determined by 

both statistical and dynamical fluctuations. Statistical 
fluctuations arise due to the finite number of particles 
measured, and can be readily calculated based on expec- 
tation values for Poisson distributions as follows: 
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The statistical fluctuations depend on the experimental 
efficiency and analysis cuts used in the reconstruction 
of charged particle trajectories (tracks). The intrinsic 
or dynamical fluctuations are deflned and evaluated as 
the difference between the measured fluctuations and the 
statistical limit 



— ,dyn — — — — ^stat- 



(3) 



As shown in [3, the dynamical fluctuations, ^dym 

can be expressed as follows: 



— Ayn — 



2R^ 
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where Rab with a, 6 = +, — are the averages of the corre- 
lation functions often used in multi-particle production 
analysis jiSlfl: 
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where i?2,afc = P2(??a, '76)/(pi,a(?7a)pi,&('7b)) - 1; Pliv) = 
dn/dr], and p2iVa,Vb) = d^n/ drjadrji, are single- and two- 
particle pseudorapidity densities respectively. The inte- 
grals could most generally be taken over the full particle 
phase space {d^p) but are here restricted (without loss of 
generality) to pseudorapidity integrals to simplify the no- 
tation. In cases where the produced particles are totally 
uncorrelated, two-particle densities can be factorized as 
products of two single-particle densities. The correla- 
tors Rab shall then vanish, and the measured dynami- 
cal fluctuations, ^dym should be identically zero. A 

deviation from zero thus should indicate correlations in 
particle production. If correlations are due to production 
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via many sub-collisions, localized sources, or clusters, one 
should further expect the strength of the correlation to 
be finite but increasingly diluted with increased number 
of production clusters or sub-collisions (hereafter called 
"clusters"). The correlators, Rat, will be inversely pro- 
portional to the multiplicity of clusters, and thus also 
inversely proportional to the total measured multiplicity 
of (charged) particles . Measurements at the ISR and 
FNAL, have shown that charged particles have long range 
(differential) correlations dominated by a dependence on 
the relative rapidity of the correlated particles. One thus 
expects, as shown in 0, that the functions, Rab, and 

^dyn should vary slowly with the detector acceptance 

as long as the rapidity width of the acceptance is smaller 
or of the order of the long range correlation width. This 
should however be experimentally verified by varying the 
acceptance used in the determination of i/^ ^dyn- 

Authors 0,01 have suggested that if the reaction dy- 
namics do not change with collision centrality, the mea- 
sure $ w {Nch)'^dyn/8 (where Nch is the charged parti- 
cle multiplicity in the rapidity range considered) should 
be independent of the collision centrality. Conversely, 
a significant collision centrality dependence of <& or re- 
lated observables should hint at a change in the collision 
dynamics. We shall thus study the collision centrality 
dependence of both ^dyii and {Nchjv^ ^dyn- The cor- 
relators, Rab, and u-^ .dym are robust variables: their 

measurements are independent of the average (global) de- 
tection efficiencies involved in the determination of mul- 
tiphcities N+ and iV_ Q. The measurement of — ^dyn 
thus does not require explicit efficiency corrections. Sec- 
ond order corrections are, in principle, needed to account 
for variations of the detection efficiency through the fidu- 
cial acceptance. In the present study, we verified that the 
relative variation of the detection efficiency (about 10% 
in the transverse momentum region under study) results 
in a systematic uncertainty less than or equal to the sta- 
tistical error of the measured values. 

The data presented are from minimum-bias and cen- 
tral trigger samples of Au -I- Au at ^Snn = 130 GeV 
acquired by the STAR experiment during the first op- 
eration of the Relativistic Heavy Ion Collider (sum- 
mer 2000). Detailed descriptions of the experiment and 
the Time-Projection-Chamber (TPC) can be found else- 
where 13]. In minimum bias mode, events were triggered 
by a coincidence between the two Zero Degree Calorime- 
ters (ZDCS) located +/- 18 m from the interaction cen- 
ter and a minimum signal in the Central Trigger Barrel 
(CTB), which consists of scintillator slats surrounding 
the TPC. The central trigger sample was acquired by 
requiring a higher multiplicity cut with the CTB corre- 
sponding to 15% of the total hadronic cross section. 

In order to minimize the need for corrections to ac- 
count for dependence of the detector acceptance and re- 
construction efficiency on the vertex position, the analy- 
sis reported here was restricted to events produced within 
±0.70 m of the center of the STAR TPC along the beam 
axis. In this range, the vertex finding efficiency is 100% 



for collisions which result in charged particle multiplici- 
ties larger than 50 tracks in the TPC acceptance. It de- 
creases to 60% for events with fewer than 5 tracks from 
the primary vertex. We verified that the measurement of 

Vj^ ^dyn is insensitive to the vertex position by comparing 

values measured for different vertex cut ranges. About 
180000 minimum-bias, and 80000 central trigger events 
were used in this analysis after cuts. 

The centrality of the collisions is estimated from the 
total charged particle track multiplicity detected within 
the TPC in the pseudorapidity range \'q\ < 0.75. We use 
eight contiguous centrality bins based on the fraction of 
triggered events: 6%, 11%, 18%, 26%, 34%, 45%, 58%, 
and 84%. The trigger efficiency is estimated to 94 ± 
2%. The above fractions thus correspond to a constant 
increase in the fraction of the geometrical cross section 
which is sampled by each multiplicity bin. 

Particle production is studied for both negative and 
positive hadrons over a transverse momentum range ex- 
tending from 0.1 to 5 GeV/c, and for pseudorapidity 
ranges from |?7| < 0.1 to 1.0 in steps of 0.1 unit of pseu- 
dorapidity. Good track quality is required by restrict- 
ing the analysis to charge particle tracks producing more 
than 15 hits within the TPC. One additionally requires 
that more than 50% of the hits be included in the final 
fit of the track. 

One uses the particle energy loss (dE/dx) measured 
with the TPC to identify the particles as pions, kaons, 
and protons (and their anti-particles). Particle identi- 
fication proceeds on the basis of a parametrization of 
the mean {{Eioss}) and width (cr) of the average energy 
loss expected for electrons, pions, kaons, and protons 
as a function of their momentum. The analyses for pi- 
ons, kaons, and protons are performed using momentum 
ranges 0.1 < p < 0.6 , 0.1 < p < 0.6, and 0.12 < p < 0.7 
GeV/c respectively. Lower bounds are set near or be- 
low detection threshold to maximize particle yields. Up- 
per bounds are used to minimize cross species contam- 
ination. The inclusive analysis of all charged species is 
performed within the range 0.1 < p < 5.0 GeV/c. Limit- 
ing the particle momenta for this analysis to less than 5 
GeV/c insured that particle charge was not mis-assigned 
while allowing for a fully inclusive measurement of the 
soft particle spectra. Given that the bulk of the parti- 
cle production is below 2 GeV/c, the inclusive analysis is 
rather insensitive to the exact value of the upper bound 
which is used. The detection efficiency rises from zero 
to roughly 85% within an interval of 0.1 GeV/c above 
detection thresholds, remaining constant for larger mo- 
menta. Measured particles are tagged as pions if their 
measured energy loss deviates by less than two standard 
deviations {2a) from the expected mean for pions of the 
same momentum, while deviating by more than 2a for 
kaons of that same momentum. Similarly particles are 
identified as kaons (protons) if the deviation from the 
kaon (proton) mean energy is less than 2a while being 
larger than 2a from the pion and proton (kaon) mean 
energy loss. Contamination of the kaons and protons by 
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pions is negligible at low momentum, and estimated to be 
less than 5% at the highest momenta accepted for those 
particles. For cross-species contamination at this level, 
it was verified that the measurement is insensitive to the 
actual value of the momentum cuts. 

To reduce contamination from secondary electron 
tracks, and focus this analysis on primary tracks, i.e. 
particles produced at the Au+Au collision vertex, only 
tracks which passed within 3 cm of the collision vertex 
were accepted. We verified electron (positron) contam- 
ination has a negligible impact on our measurements of 

_dyn by repeating the analysis with and without an 

electron/positron exclusion cut based on the track energy 
loss measured in the TPC, i.e. accepting tracks with a 
dE/dx more than 2 standard deviations away from the 
expected value for an electron of the measured momen- 
tum. 

As already mentioned, the measurement of i^^ _dyn is 

independent of the average detector efficiency. It is there- 
fore also insensitive to particle losses, e.g. anti-protons, 
due to scattering through the detector. It is however 
sensitive, in principle, to the generation of background 
particles within the detector. The effect of such back- 
ground particles (e.g. protons scattered off the beam 
pipe) is minimized by using the 3 cm distance of closest 
approach cut mentioned above. Also, it was considered 
whether finite track splitting, possibly encountered in the 
reconstruction of charged particle tracks in the TPC, may 

produce measurable effects on ^dyn- We verified that, 

within statistical errors, the same value is obtained when 
the pseudorapidity regions used to count positive and 
negative tracks were separated by a Arj = 0.25 gap. 

Since finite width multiplicity bins were used for this 

analysis, values of i/^ ^dyii are multiplicity-bin averaged 

according to the following expression: 



(6) 

where P{M) is the probability of having a total charge 
multiplicity, M, and ,dyn{M) is given by 



V+^^dyniM) = 



{N+{N+ - 1))m 
{N-{N- - 1))m 



(7) 



IM 



{N+)m{N-)m 

The notation < O >m is used to indicate the average of 
the quantity O for all events with a charged particle mul- 
tiplicity, M, in the pseudorapidity range \ri\ < 0.75. Our 
analysis proceeds in two passes. The first pass involves 
the determination of the averages {N±)m as a function 
of the multiplicity, M, using unity bin- width in M while 
the second pass uses these averages as coefficients in the 
above expression of h'^^-^^dyniM). The averages {N±)m 



are determined from the events with multiplicity M: 

The sum is taken over the Nev{M) events of multiplic- 
ity, M, present in our sample. The averages {N±)m thus 
obtained display a scatter determined by the finite statis- 
tics about a monotonically increasing trend (with M). If 
uncorrected, this scatter, may induce an artificial change 

of the value of i'^ ,dyn(M) in each bin. To minimize this 

effect, we model (fit) the average {N±)m dependence on 
the multiplicity M with a polynomial optimized to repro- 
duce the shape of the dependence. We then determine 
iy+-^dyn{M) using the averages {N±)fit.M = N±m pre- 
dicted by the fit rather than the actual averages. The 

calculation of u^^ _dyn in a- finite width multiplicity bin 

then proceeds with the following expression: 



V+~,dyn (A/i„,„<Af<M,„gh) 



(9) 
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N+{N+ - 1) 



N\ 



+ .M 



N_{N_ - 1) 
-2_^^^- 



.M 



where the sum is taken over the Ng^ events in the multi- 
plicity bin Mion, < M < Muigh- 

The quantity {N)v^ ^dyn is determined in a similar 

fashion using the following expression: 



{N)v+^^dyn (Mia„<M<Mh,gh) = 



(10) 



+ 



events 

N+{N+ - 1) 
iV_(A^_ - 1) 



.M 



To study the effect of this method of bin-averaging, a sim- 
ulation was performed using HIJING events, comparing 
the results of equation 10, and equation 3 in the limit of 
large statistics. The HIJING model does not incorporate 
re-scattering and should not therefore exhibit a signifi- 
cant centrality dependence. The results showed that for 
all bins except the lowest multiplicity bin used for this 
analysis, the two equations gave the same result within 
the quoted systematics. In the first multiplicity bin, 
equation 10 yielded a result approximately 15% larger 
than equation 3. 

Fig. ^ shows the dynamical fluctuations, ,dyn, 

of the net charge measured in the pseudorapidity range 
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\ri\ < 0.5, as a function of the total multiplicity, M, mea- 
sured in the pseudorapidity range \ri\ < 0.75. The hor- 
izontal bars on the data points reflect the width of the 
multiplicity bins used in this analysis while the vertical 
bars reflect statistical errors. We estimate the systematic 
errors based on data taken and analyzed with different 
trigger and analysis cuts, to be of the order of 2%. An 
additional systematic uncertainty of the order of 3% is 
derived by a separate analysis of different data subsets. 
The dynamical fluctuations of the 5% most central colli- 
sions then amount to ^dyn = — 0.00236±0.00006(stat) 

±0.00012(syst). The dynamical fluctuations are finite 
and negative: a clear indication that positive and nega- 
tive particle production are correlated within the pseudo- 
rapidity range considered (see Eq. One observes the 
strength of the dynamical fluctuations decreases mono- 
tonically with increasing collision centrality. This can be 
understood from the fact that more central Au-|-Au col- 
lisions involve an increasing number of "sub-collisions" 
(e.g. nucleon-nucleon collisions): the two-particle corre- 
lations are thus increasingly diluted and the magnitude 

of ^dyn is effectively reduced. 

We compare our results, for the most central colli- 
sions, to those recently reported by the PHENIX col- 
laboration |l4j which measured net charge fluctuations 
in terms of the relative variance ujq = {AQ'^) /Nch in 
the rapidity range \ri\ < 0.35, and the angular range 
A$ = 7r/2, for p±_ > 200MeV/c. They reported a value 
iVQ = 0.965 ± 0.007(stat) -0.019(syst) for the 10% most 
central collisions. The large (unidirectional) systematic 
error is reported to result from uncertainty in the correc- 
tion applied for effects of flnite detector efficiencies. In 
order to compare the PHENIX result with the present 
study, we use the expression in reference 
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The charged particle multiplicity in the PHENIX detec- 
tor acceptance is 79 ± 5 for the 10% most central col- 
lisions. This comparison gives i^^ _dyn — —0.0018 ± 

0.0004 {stat) — 0.009 {syst) in agreement with the value 
of v+-,dyn = -0.00263 ± 0.00009(stat) ±0.00012(syst) 
we measure for 11% central collisions. The agreement 
is best if one considers the low bound of the PHENIX 
measurement which is maximally corrected for finite ef- 
ficiency (which is refiected in the systematic error). The 
difference between the two results might be due, in part, 
to dependence of the multiplicity fiuctuations on rapidity 
and azimuthal angle as well as acceptance effects. 

It is important to consider the effects of charge con- 
servation on the net charge fluctuations since they are 
expected to be non-negligible even for small finite rapid- 
ity coverage 0. The contribution is estimated to be 
— 4/(iV)4^ where (iV)4^ is the total number of charged 
particles produced by the collisions. The PHOBOS col- 
laboration has reported [ll| that the total charged parti- 
cle multiplicity amounts to 4200±470 in the 6% most cen- 
tral Au-|-Au collisions at y^sj^ =130 GeV. The charge 



FIG. 1: (a) Dynamical fluctuations, ,dyn, measured in 

\ri\ < 0.5 as a function of the collision centrality estimated 
with the total (uncorrected) multiplicity, M, in \ri\ < 0.75. 

(b) {N)v^ ^dyn measured in \ri\ < 0.5 vs M (opened circles) 

compared to the charge conservation limit (dotted line), res- 
onance gas expectation based on ref. 0] (solid line); and 
HIJING calculation (solid squares) 



conservation contribution to the measured dynamical 
fiuctuations is thus of the order of -0.00095 ± 0.0001, 
i.e. 40% of the observed dynamical fluctuations. 

We next discuss the centrality dependence of the fluc- 
tuations. In central collisions, the measured dynamical 

fluctuations, ^dyn are expected to be reduced due to 

dilution of the two-particle correlations. One expects 

the magnitude of v-^ ^dyn should scale inversely to the 

number of sub-collisions producing particles. Assuming 
the average number of particles produced by such sub- 
collisions is independent of the collision centrality, one 
then expects the fluctuations to scale inversely as the 

charged particle multiplicity. The quantity {N)v^ ^dyn 

should therefore be independent of collision centrality if 
no significant variation in the mechanism of the parti- 
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FIG. 2: Fluctuations v^-^dyn for the 6% most central colli- 
sions as a function of the range of integrated pseudorapidities. 
The expected limit due to charge conservation is shown as a 
dotted line. 



cle production arises with collision centrality. This no- 
tion was suggested by Gazdzicki j3| and Mrowczyn- 
ski in terms of the fluctuation measure <i> which, as 
shown in 0, is equal to {N)v+^^dynl^ for {N+) w (A^-). 
Fig. shows the measured centrality dependence of 
(7V)i/_| ,iym calculated with Eq. ^| for all charged par- 
ticles produced in the pseudorapidity range jryj < 0.5. In 
this Figure, the charged particle multiplicity, N ^ is cor- 
rected for finite detection efficiencies using correction fac- 
tors which depend linearly on the charged particle mul- 
tiplicity (TPC detector occupancy) with values ranging 
from 85% to 70% for peripheral and central collisions re- 
spectively [T^ . The measured values range from —1 to 
— 1.4 and are approximately a factor of two larger than 
the charge conservation limit, shown as a dotted line, in 
Fig. ^p). This indicates dynamical fluctuations are not 
only finite but in fact rather large. As discussed in de- 
tail below, the values measured for {N)v^ ^dyn however 

fall short of predictions for a resonance gas in equilib- 
rium (approximately —1.7; solid line) and for a scenario 
involving a quark-gluon gas undergoing fast hadroniza- 
tion (approximately -3.5; not shown in Fig. QJ?) 
The measured values are in qualitative agreement with 
a calculation based on HIJING (solid squares) I^i" 
deed, the values predicted by HIJING are within 20% 
of the measured values at all centralities. While the HI- 
JING calculation is independent of collision centrality, 
the experimental data exhibit a small but finite central- 
ity dependence which is significant above the first bin 
in Fig. lb. The HIJING calculation does not feature 
rescattering, and is therefore not expected to exhibit a 
significant centrality dependence. The observed central- 
ity dependence may then suggest there are rescattering 
effects, or other dynamical effects with centrality, and its 
interpretation requires further investigation. 

The magnitude of the net charge dynamical fluctua- 
tions is determined by the strength of the two-particle 
correlations in the integrated rapidity range. Measure- 



ments from p-f p collisions at the ISR and p + p collisions 
at FNAL indicate that the relevant rapidity interval for 
two-particle correlations is approximately one unit. One 
thus expects the dynamical fluctuations to exhibit a mild 
dependence on the rapidity range used for the measure- 
ment 7]. Fig. 121 shows the measured dynamical fluctu- 
ations (filled circles) as a function of the pseudorapidity 
range. The pseudorapidity integration range is varied 
from —0.1 < 7y < 0.1 to —1.0 < rj < 1.0 in discrete steps 
of 0.1 units of pseudorapidity. Error bars shown are sta- 
tistical only. Focusing on the region in Fig. 2 where 
systematic effects due to finite are expected to be small, 
we examine the data for rj > 0.4. One observes the abso- 
lute value of the the dynamical fluctuations is largest in 
this range for |?7| 0.4, and that it decreases monoton- 
ically for larger acceptance without, however, reaching 
the charge conservation limit. One finds |t'dyn| decreases 
by 35-40% while the integrated pseudorapidity range is 
increased by a factor of 5 from 0.4 to 2 pseudorapidity 
units. The dependence of dynamical fluctuations on the 
experimental acceptance is rather modest. In contrast, 
the $ measure increases approximately by a factor of 10 
from —0.1 < 1] < 0.1 to —1.0 < rj < 1.0 due to its explicit 
dependence on the pseudorapidity bin size. 

We next consider the above results in the light of cor- 
relation functions measured in p+p and p+p collisions at 
CERN and FNAL H El [H with the use of Eq. H To 
account for the unavailability of p-l-p comparison data at 
the same energy as RHIC, an interpolation was made us- 
ing results obtained at lower and higher collision energies 
(parameterization from jl^l). Based on results published 
in Refs H, 0, ^ , we also note that the correlation func- 
tion for oppositely charged particles, (y+ « ?/_), is 

found to be approximately twice as strong as the same 

sign particles correlations, ~ R [1,0], and that 

it is independent of the collision energy. The CERN and 
FNAL measurements H 0, find the single charged 
particle and two-particle (charged-charged) pseudorapid- 
ity densities to be respectively pi{r] = 0) « 2.06 and 

^2(0,0) = P2{m = o,?72 = 0) - pi[m = 0)^1(772 - 

0) « 2.8. The charged-charged correlation integral 
Rcc = {R++ + R-- + 2i?+-)/4 is thus Rcc ~ 0.66 (see 
ref. 0). Furthermore, assuming equal multiplicities of 
positively and negatively charged particles, one finds for 
the charged-charged correlation Rcc ~ l■5R-^.^, which we 
use to estimate the correlation measured in this work as 
R++ + R— - 2i?+_ « -2R++ « 4i?cc/3 « 0.88. The 
pseudorapidity densities are very different in p-l-p and 
A+A collisions. Under assumption that the correlations 
are due to production in a finite number of sources (clus- 
ters) , they should be inversely proportional to the parti- 
cle density. In the 5% most central Au-|-Au collisions, the 
pseudorapidity charged particle density {dN/ drf) is about 
526 ± 2{stat) ± 36{syst) [l^ compared to approximately 
2.06 in pp collisions. Such a dilution would give for the 
correlation function a value of 0.88 • 2.06/526 « 0.0034, 
in qualitative agreement with the measured values for 
Au-|-Au collisions presented in this paper. We stress that 
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valuable insight can be gained by comparing the current 
130 GeV data and upcoming 200 GeV Au+Au analysis 
with explicit measurements made in p+p collisions rather 
than using the above first order approximation. 

We next compare our measurement of the dynami- 
cal fluctuations to predictions of net charge fluctuations 
based on thermal models H H IH HI H^. To this 

end, we express our measurement of ^dyn in the range 

\r]\ < 0.5 in terms of the D variable introduced in [3, 
using 

D = 4+{N)iy+^^dvn (12) 

valid for N+ ^ N_ 0. We find using data shown in 
Fig. ^ that D decreases from 3.1 ± 0.05 (statistical er- 
ror only) for the most peripheral collisions measured to 
2.8 ± 0.05 in central collisions. However, a comparison 
to thermal model predictions requires the data to be cor- 
rected for charge conservation effects. One must subtract 
the charge conservation contribution which amounts to 
AD = -0.00095 X 526 = -0.50±0.06. The corrected val- 
ues of D thus range from 3.6±0.1to3.2±0.1. According 
to the discussion of Refs. H S IH IH 113 , these values 
approach that {D « 2.8) expected for a resonance gas. 
They are significantly l arg er than expected in the above 
referenced work '5','2l',|2l|2l| for a quark-gluon gas un- 
dergoing fast hadronization and freeze-out (Z? « 1). It is 
not possible to draw a firm conclusion concerning the ex- 
istence or non-existence of a deconfined phase during the 
collisions from these results since, as the above authors 
have pointed out, incomplete thermalization could lead 
to larger fluctuations than expected for a QGP. Other 
work |2J| has also suggested that the prediction of I? w 1 
for a quark-gluon gas is model dependent, and that other 
effects such as gluon fragmentation prior to to hadroniza- 
tion could increase the fluctuations expected even if a 
quark-gluon plasma were produced. 

We extend the study of net charge fluctuations to iden- 
tified particles and consider measurements of the net 
charge fluctuations of pions, kaons, and protons/anti- 
protons. Measurement of the K'^,K~ andp,p net charge 
are of particular interest as they address respectively fluc- 
tuations of net strangeness and baryon number which 
might be more sensitive to the details of the collision pro- 
cess. The results are compiled in Table 1 for all charged 
species, pions, kaons, and p,p. The results indicate 
that the dynamical fluctuations for pions are approxi- 
mately the same magnitude as for inclusive non-identified 
charged particles. One however discerns a small but finite 
difference, especially for integrated pseudorapidity ranges 
\ri\ < 0.7 and larger. The measurement of K^,K~ and 
p,p fluctuations is hampered by the smaller multiplici- 
ties and finite detection efficiencies for kaons and pro- 
tons and their anti-particles. Our measurement, which 
is presented in Table 1 for acceptances from I77I < 0.5 
to 1 77 1 < 1.0 is thus limited to a central collision trig- 
ger sample. The effect of the variation of efficiency near 
detection threshold was studied by changing the trans- 
verse momentum threshold used in the determination of 



^ciyn- It was found that for inclusive non-identified 

particles, ,dyn changed by less than 3% while varying 

the transverse momentum cutoff for particle detection 
from 0.1 to 0.2 GeV/c. The same study using HIJING 
events led to a 10% change in ^dyn- 

The systematic error for protons (anti-protons) is dif- 
ficult to assess, since GEANT studies indicate a consid- 
erable fraction of the proton yield below 0.4 GeV/c is 
associated with pion-induced proton knockout reactions 
in the beam pipe. Background protons bear little cor- 
relation with anti-protons. The terms R++ and 

involved in the calculation of ,dyn should have a Pois- 

sonian behavior, and therefore the contribution of uncor- 
related background to these terms should partly cancel. 

We find the value of i'^ ^dyn exhibit changes smaller than 

the statistical uncertainties when raising the threshold 
from 0.2 to 0.3 GeV/c, and hence ascribe a systematic 
error of the order of 20% for the p,p measurement. 

The dynamical fluctuations of the charged kaons and 
p,p are also finite. Their size (absolute value) are in 
fact larger than the dynamical fluctuations measured for 
pions and for inclusive non-identified charged particles. 
The proton dynamical fiuctuations are somewhat larger 
than the kaon fluctuations. Strangeness conservation, 
and baryon number conservation should influence the size 
of the dynamical fluctuations for the net charge of kaons, 
and p,p respectively. The charge conservation limit de- 
rived for inclusive non-identifled charged particles can be 
readily reinterpreted to estimate the expected magnitude 
of dynamical fluctuations for K~^,K~ and p,p. One flnds 
that the kaon and pjj dynamical fluctuations are similar 
or slightly larger than their respective charge conserva- 
tion limits. 

TABLE I; 1000 x v+-,dyn, for charge pions, kaons, and pro- 
tons, as a function of the integrated pseudorapidity range 



H 


All +- 


7r± 




P,P 


0.5 


-2.36 ±0.06 


-2.4 ±0.1 


-5 ±3 


-3 ±7 


0.6 


-2.27 ±0.04 


-2.4 ±0.1 


-5 ±2 


-5 ±3 


0.7 


-2.11 ±0.04 


-2.18 ±0.08 


-4 ±2 


-7 ±5 


0.8 


-1.98 ±0.03 


-2.12 ±0.07 


-6 ±2 


-8 ±3 


0.9 


-1.90 ±0.03 


-2.02 ±0.06 


-6 ±2 


-9 ±2 


1.0 


-1.75 ±0.02 


-1.92 ±0.06 


-7± 1 


-8 ±2 



Summary and Conclusions 

We have measured event-by-cvent net charge dynami- 
cal fluctuations for inclusive non-identifled charged par- 
ticles, as well as for identifled pions, kaons, and protons 
and their anti-particles in Au±Au collisions at ^/snn ~ 
130 GeV. Dynamical fluctuations measured for inclusive 
non-identifled charged particles are flnite and exceed by 
nearly a factor of two expectations based on charge con- 
servation. We flnd the magnitude of the net charge dy- 
namical fluctuations to be in qualitative agreement with 
expectations based on measurements of charged parti- 
cle correlation functions in p+p collisions measured at 
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the ISR. We however find that although the fluctuations 
roughly scale in proportion to the reciprocal of the pro- 
duced charged particle multiplicity, the scaling is not per- 
fect, and the quantity {N)i>^ iiyn exhibits a small de- 
pendence on collision centrality, which suggests the two- 
particle correlations may be modified in central collisions 
relative to peripheral collisions. 

A comparison of our measurement with thermal model 
predictions 0, |^ appear to indicate fluctuations 
at a level that might be expected if the Au-|-Au system 
behaved like a resonance gas. Although the size of the 
fluctuations is significantly larger than expected in that 
work for a quark-gluon gas, limitations of the model used 
prevent a conclusion on the existence or non-existence of 
a quark-gluon plasma phase based on these results. 

Finally, we report the first measurement of net-charge 
dynamical fiuctuations of identified pions, kaons, and 
protons. Pions exhibit dynamical fluctuations slightly 
larger than the values obtained with our inclusive mea- 
surement. Kaons and protons are found to exhibit dy- 
namical fluctuations that are 2 to 4 times larger than 
those observed for all charged particles. However, the 
lower production multiplicities of these particles may im- 
ply the dynamical fluctuations are dominated by charge 
conservation effects. Further data are needed to assess 
whether the dynamical fluctuations of kaons (protons) 



significantly exceed the minimal values constrained by 
strangeness (baryon) charge conservation. 
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